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Abstract 

We investigate Bose-Einstein condensation of noninteracting gases in a harmonic trap with an 
off-center dimple potential. We specifically consider the case of a tight and deep dimple potential 
which is modelled by a point interaction. This point interaction is represented by a Dirac delta 
function. The atomic density, chemical potential, critical temperature and condensate fraction, 
the role of the relative depth and the position of the dimple potential are analyzed by performing 
numerical calculations. 

Topic: Physics of Cold Trapped Atoms 

Report number: 6.13.1 



Electronic address: dtarhan@gmail.com 



I. INTRODUCTION 



The phase space density of a Bose-Einstein condensate (BEC) can be increased by mod- 
ification of the shape of the potential jl. "Dimple" -type potentials are the most commonly 
used potentials for this purpose [2I, Is,!^. A small dimple potential at the equilibrium point 
of the harmonic trapping potential is used to enhance phase-space density by an arbitrary 
'actor ^. A tight dimple potential is used for a recent demonstration of caesium BEC a 
sl- Quite recently, such potentials were proposed for efficient loading and fast evaporative 
cooling to produce large BECs [sl- Attractive applications, such as controlling interaction 
between dark solitons and sound [6|, introducing defects such as atomic quantum dots in 
optical lattices [?], or quantum tweezers for atoms [Is] are offered by using tight dimple poten- 
tials for (quasi) one-dimensional BECs . Such systems can also be used for spatially selective 
loading of optical lattices 9|. In combination with the condensates on atom chips, tight and 
deep dimple potentials can lead to rich novel dynamics for potential applications in atom 
lasers, atom interferometers and in quantum computations (see Ref. ^Oi] and references 
therein) . 

In this paper we continue the discussion of our recent paper In that paper, we 

modelled the dimple type potentials by Dirac 6 functions and investigated the change of 
chemical potential, critical temperature and condensate fraction of a harmonic trap with 
respect to the various strengths of Dirac 6 functions. In this paper, we investigate the 
behavior of the same physical quantities for a 6 function which can be located at different 
positions than the center of a harmonic trap. We find that while a centrally positioned 
dimple potential is most effective in large condensate formation at enhanced temperatures, 
there is a critical location for which the condensate fraction and the critical temperature 
can also be enhanced relatively. This might be useful in spatial fragmentation of atomic 
condensates. 

The paper is organized as follows. In Sect. II, we review shortly the analytical solutions of 
the Schrodinger equation for a harmonic potential with a finite number of Dirac 5-decorated 
harmonic potential and give the eigenvalue equation of the harmonic potential with a Dirac 
S function. In Sect. Ill, determining the eigenvalues numerically, we show the effect of the 
dimple potential on the condensate fraction and the transition temperature and investigate 
the change of this values with respect to the position of the Dirac 6 function. Finally, we 
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conclude in Sect. IV. 



II. HARMONIC POTENTIAL DECORATED WITH DIRAC DELTA FUNC- 
TIONS 



We begin our discussion by reviewing the one dimensional harmonic potential decorated 
with the Dirac 6 functions lll-jl^ for the seek of completeness. The potential for is given 



as: 



1 ^ 

V{x) = -mu^x^ ~" 2m ^ ~ 



(1) 



where u is the frequency of the harmonic trap, P is a finite integer and cTj's are the strengths 
(depths) of the dimple potentials located at Xj's with xi < X2 < ... < xp with Xi G (— oo, oo). 
The factor /2m is used for calculational convenience. Negative (Ti value represents repulsive 
interaction while positive (Ti value represents attractive interaction. The time-independent 
Schrodinger equation equation for this potential is written as: 



+ V{x)^!{x) = E^{x). 



(2) 



2m dx'^ 

By defining E = {C^ + j)huj, with ^ a real number, and introducing dimensionless quantities 
z = x/xq, and Zi = Xi/xq with xq = ^Jfi/2muo, the natural length scale of the harmonic 
trap, we can re-express Eq. ([2]) as 

9 P 



d?m{z) 
dz^ 



+ 



^Ai5{z - Zi) 

i 

n 



2 4 

where Aj = Xocxj. By using transfer matrix approach 
eigenvalue equation: 



^(z) =0, 



(3) 



16| . we get the following 



^ X,D^{z^) D^{-z,) _ ^ 
W 



(4) 



for ^ and using E = + ^^I^jJ- Here, D^{z) and D^{—z) are parabolic cylinder functions 
and Zi = Xi/xq. The Wronskian W of D^{z) and D^{—z) is 

2(«+3/2)^ 



W = W[D^{z),D^{-z)] 



For zi = 0, these results reduce to the results in Ref. 



r(-|)r(¥) 

iil,Q,Q. 



(5) 
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III. BEC IN A ONE-DIMENSIONAL HARMONIC POTENTIAL WITH A DIRAC 
S FUNCTION 



In this section we calculate the condensate fraction, chemical potential, critical tem- 
perature and density profile for different depths, sizes and positions of a dimple potential 
modelled by a Dirac 6 function. In order to describe the depth and size of a dimple potential 
in a systematic way we define a dimensionless variable in terms of the strength of the Dirac 
6 functions as: 

We will present our results with respect to A and zi defined in the previous section. 



In ref . 



Ill], we have estimated a values approximately according the parameters in refs. 



17j and [sj. We find that, if 10^1/m < a < lO^^l/m then 320 < A < 32000 for the 
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and for the experimental 



we have shown that, 



experimental parameters m= 23 amu (^^Na), = 27r x 21 Hz 
parameters m= 133 amu (^'^^Cs), u; = 27r x 14 Hz [sl. In ref. 
even for small A values, condensate fraction and critical temperature change considerably. 
How sensitive such a change would occur depending on the location of the dimple trap was 



a question left unanswered in Ref. 



ll|. 



We begin our discussion by investigating the change of the critical temperature with 
respect to the position of Dirac 6 function. The critical temperature (Tc) is obtained by 
taking the chemical potential equal to the ground state energy [n = Eg = Eq) and by 
solving 

i=l 

for Pc where Pc = ^/{kBTc). For finite value, we define T° as the solution of Eq. ([7]) for 
A = (only the harmonic trap). 

In Eq. (171), ej's are the eigenvalues for the harmonic potential decorated with a single 
eccentric dimple potential at zi. The energies of the decorated states are found by solving 
Eq. (jl]) numerically. Then, these values are substituted into the Eq. ([7]); and finally this 
equation is solved numerically to find T^. We obtain for different Zi values and A = 32 
and show our results in Fig. ([1]). Since harmonic potential is symmetric negative zi values 
will give a the same values for critical temperature with positive ones. As zi increases, 
the energy of the ground state increases so that the critical temperature decreases as zi 
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FIG. 1: The critical temperature Tc for A'^ = 10^. A is a dimensionless variable defined in Eq. 
Here we use m = 23 amu (^^Na) and w = 27r x 21 Hz {l?! . 



gets larger. On the other hand, as the dimple trap becomes farther to the center of the 
harmonic trap, the critical temperature cease to decrease and starts rising again as seen in 
Fig. ([1]) around Zi = 3.5. Finally, at very large separations between the dimple trap and the 
harmonic trap center, the critical temperature no longer changes with the location of the 
dimple trap and saturates at the value corresponding to the of the critical temperature for 
the single harmonic trap per se. The increase of the critical temperature around zi = 3.5 can 
be explained as follows: As Zi increases it becomes closer to the node of the first exited state 
wave function of the harmonic potential. At that value the change of the energy eigenvalue of 
the first excited state vanishes and the difference between the first excited state and ground 
state eigenvalues increase. Thus, the particles favor the ground state which increases the 
critical temperature. In Fig. ([T]) we take N = 10^ and use typical experimental parameters 
m= 23 amu (^^Na) and = 27r x 21 Hz [17|. 

For a gas of N identical bosons, the chemical potential /i is obtained by solving 



E 



E 



1 



(8) 



i=0 i=l 

at constant temperature and for given N, where is the energy of state i. We present the 
change of as a function of T/T° in Fig. (El) for A^ = 10^; A = 0; A = 32, = and A = 32 
Zi = 1. 
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FIG. 2: The chemical potential vs temperature T/T^ for N = 10"^ . The solid line shows fj, only 
for the harmonic trap. The dotted line shows /x for zi = 1 and A = 32. The dashed line shows /U 
for zi = and A = 32 . The other parameters are the same as Fig. [TJ 

By inserting /i values into the equation 

we find the average number of particle in the ground state. Nq/N versus T /T^ for = 10^ 
A = 32 are shown in Fig. ([3]). In this figure, the solid line shows the condensate fraction for 
zi = and the dashed line shows the condensate fraction for zi = 1. In Ref. 18|, Ketterle et. 
al mentions that the phase transitions due to discontinuity in an observable macro parameter 
occurs only in thermodynamic limit, where N oo. However, we make our calculations for 
a realistic system with a finite number of particles in a confining potential. Thus, N^/N is 
a finite non-zero quantity for T < Tc without having any discontinuity at T = Tc. 

We also investigate the behavior of the condensate fraction as a function of the position 
of Dirac 6 for A = 32, T = and present the results in Fig. (jlj). 

Finally, we compare density profiles of condensates for a harmonic trap and a harmonic 
trap decorated with a delta function (A = 3.2 and zi = 1) in Fig. Since the ground state 
wave functions can be calculated analytically for both cases, we find the density profiles by 
taking the absolute square of the ground state wave functions. Comparing the graphics of 
density profiles, we see that an offcenter dimple potential maintain a higher density at the 
position of the Dirac 6 function which may be utilized for the fragmentation of a BEC. 



u 

1 2 

c 



FIG. 3: Nq/N vs T/T^ for = 10^ and A = 32. The solid line for zi = 0, the dashed line for 
zi = I. The other parameters are the same as Fig. [H 




FIG. 4: Nq/N vs r/r° for N = 10^ and A = 32. The other parameters are the same as Fig. [H 



IV. CONCLUSION 



We have investigated the effect of the location of the tight dimple potential on the results 
reported recently in our paper ll|]. We model the tight dimple potential with the Dirac 
6 function. This allows for analytical expressions for the eigenf unctions of the system and 
simple eigenvalue equations greatly simplifies subsequent numerical treatment. We have 
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FIG. 5: Comparison of density profiles of a BEC in a harmonic trap with a BEC in a harmonic trap 
decorated with a 6 function (A = 3.6, zi = 1). The solid curve is the density profile of the BEC 
in decorated potential. The dashed curve is the density profile of the ID harmonic trap (A = 0). 
The parameter z is dimensionless length defined after Eq. ([2]). The other parameters are the same 
as Fig. [H 

calculated the critical temperature, chemical potential, condensate fraction and presented 
the effects of the location of the dimple potential. We find that the dimple type potentials 
are most effective when they are applied to the center. While it is also advantageous to 
place the dimple potential at the nodes of the excited state, where our results revealed a 
relative enhancement of the critical temperature and the condensate fraction. Determining 
the density profiles of the BECs in the harmonic trap and in the decorated trap with the 
Dirac 6 function at this critical position, we argued that eccentric dimple trap at such a 
critical location can be used for spatial fragmentation of large, enhanced BECs. 

The presented results are obtained for the case of noninteracting and one dimensional 
condensates for simplicity. In such a case, stability of the condensate may become question- 
able and should be addressed separately in detail, [lol. The treatment should be extended 
for the case of interacting condensates in larger (or quasi) dimensional traps in order to 
make the results more relevant to experimental investigations. 
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Figure Captions 

Fig. 1 The critical temperature Tc for = 10^. A is a dimensionless variable defined in 
Eq. (p. Here we use m = 23 amu (^^Na) and = 27r x 21 Hz Q- 

Fig. 2 The chemical potential vs temperature T/T° for = 10"^ . The solid line shows 
/i only for the harmonic trap. The dotted line shows fi for Zi = 1 and A = 32. The dashed 
line shows fi for zi = and A = 32 . The other parameters are the same as Fig. [H 

Fig. 3 No/N vs T/T° for = 10^ and A = 32. The solid hue for zi = 0, the dashed 
line for Zi = 1. The other parameters are the same as Fig. [H 

Fig. 4 No/N vs T/T° for = 10^ and A = 32. The other parameters are the same as 

Fig. m 

Fig. 5 Comparison of density profiles of a BEC in a harmonic trap with a BEC in a 
harmonic trap decorated with a 6 function (A = 3.6, zi = 1). The solid curve is the density 
profile of the BEC in decorated potential. The dashed curve is the density profile of the ID 
harmonic trap (A = 0). The parameter z is dimensionless length defined after Eq. ([2]). The 
other parameters are the same as Fig. [H 
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